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Abstract 

We calculate the thermodynamic functions of a quark-gluon plasma for general N c and Nf to 
three-loop order using hard-thermal-loop perturbation theory. At this order, all the ultraviolet 
divergences can be absorbed into rcnormalizations of the vacuum, the HTL mass parameters, and 
the strong coupling constant. We show that at three loops, the results for the pressure and trace 
anomaly are in very good agreement with recent lattice data down to temperatures T ~ 2 T c . 
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1. Introduction 

The ultrarelativistic heavy- ion collision experiments at Brookhaven National Labs (RHIC), and 
CERN (LHC) allow the experimental study of matter at energy densities exceeding that required 
to create a quark-gluon plasma. At RHIC, the initial temperatures were up to twice the critical 
temperature for deconfincmcnt, T c ~ 170 McV. This corresponds to a strong coupling constant of 
a s = g 2 s /4ir ~ 0.3. Theoretically, one expected that this state of matter could be described in terms 
of weakly interacting quasiparticles; however, data from RHIC suggest that the state of matter 
created behaves more like a strongly coupled fluid with a small viscosity [l[ . This has inspired work 
on strongly-coupled formalisms based on e.g. the AdS/CFT correspondence. 

In the upcoming heavy- ion collisions at LHC, the energy densities and therefore the initial 
temperatures will be higher than those at RHIC. One expects temperatures up to 4-6T c and due 
to asymptotic freedom of QCD, this corresponds to a smaller coupling constant. An important 
question is then whether the matter generated can be described in terms of weakly interacting 
quasiparticles at these higher temperatures. Lattice simulations of QCD provide a clean testing 
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ground for the quasiparticle picture and in this Letter we compare new next-to-next-to-leading order 
(NNLO) results for thermodynamic functions of QCD with lattice data [D, 0] and with previous 
results at leading order (LO) and next-to-leading order (NLO) The calculation is based on 
hard-thermal-loop perturbation theory (HTLpt) which is a reorganization of finite-temperature 
perturbation theory. In HTLpt one expands around an ideal gas of massive gluonic and quark 
quasiparticles where screening effects and Landau damping are built in. Our results indicate that 
the lattice data are consistent with a quasiparticle picture down to temperatures of T ~ 2 T c , 
depending on which thermodynamic function one considers. 

The calculation of thermodynamic functions for quantum field theories at weak coupling has a 
long history. The free energy of QCD is now known up to order a;?loga s 0, @. Unfortunately, a 
straightforward application of perturbation theory is of no quantitative use at phenomenologically 
relevant temperatures. The problem is that the weak-coupling expansion oscillates wildly and shows 
no sign of convergence unless the temperature is astronomically high. For example, if one compares 
the (^-contribution to the QCD free energy with three quark flavors to the gf-contribution, the 
former is smaller only if a s < 0.07, which corresponds to T ~ 10 5 GeV or T ~ 5 x 10 5 T c . 

There are several ways of reorganizing the perturbative series at finite temperature Q and they 
are all based on a quasiparticle picture where one is perturbing about an ideal gas of massive 
quasiparticles, rather than that of an ideal gas of massless quarks and gluons. In scalar <^> 4 -theory 
the basic idea is to add and subract a thermal mass term from the bare Lagrangian and to include 
the added piece in the free part of the Lagrangian. The subtracted piece is then treated as an 
interaction on the same footing as the quartic term Q. In gauge theories, however, simply adding 
and subtracting a local mass term, violates gauge invariance [9(- Instead, one adds and subtracts 
an HTL improvement term, which dresses the propagators and vertices self-consistently so that the 
reorganization is manifestly gauge invariant [l0( . 



2. Hard-thermal-loop perturbation theory 

The Lagrangian density for an SU(iV c ) Yang-Mills theory with Nf fermions in Minkowski space 

is 

£qcd = -^Tr[G MI ,GH + #7 M -D/^ + £gf + £ gh + A£ QCD . (!) 

where the field strength is = d^A" — d v A^ — ig^A^, A v \ and the covariant derivative is 

D 1 * = d^—ig s A^. A£qcd contains the counterterms necessary to cancel the ultraviolet divergences. 
The ghost term £ g h depends on the gauge- fixing term £ g f. In this paper we choose the class of 
covariant gauges where the gauge-fixing term is 



{d^f . (2) 



HTLpt is by construction gauge invariant order by order in perturbation theory and our results 
are therefore independent of the gauge- fixing parameter £. In Ref. the gauge- fixing parameter 
independence in general Coulomb and covariant gauges was explicitly demonstrated at NLO. Fur- 
thermore, we use MS dimensional regularization with a renormalization scale to regularize infrared 
and ultraviolet divergences. With the standard normalization, we have ca = N c , cLa = N% — 1, 
s F = Nf/2, d F = N c N f , and s 2F = (iV c 2 - l)N f /AN c . 
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Hard-thermal-loop perturbation theory is a reorganization of the perturbation series for thermal 
QCD. The Lagrangian density is written as 



C = (£q C d + £htl) _ + ACutl , (3) 

g s ^VSg s 

where A£htl contains the additional countcrterms necessary to cancel the ultraviolet divergences 
introduced by HTLpt. The HTL improvement term is 



£htl = 6)m 2 D Tr (c, a (jfj^) G% j + (1 - 5) im\^ (j^j $ , (4) 



v 

where y M = (l,y) is a light-like four-vector, and {. . .) y represents the average over the directions of 
y. The free parameters mjj and m q are identified with the Debye screening mass and the fermion 
thermal mass. The parameter 5 is a formal expansion parameter and bookkeeping device: HTLpt 
is defined as an expansion in powers of S around S = 0. This expansion generates systematically 
dressed propagators and vertices. It also automatically generates new higher-order terms that ensure 
that there is no overcounting of Feynman diagrams. The HTL perturbative expansion generates 
ultraviolet divergences. There is no general proof that HTLpt is renormalizable, so the general 
structure of the counterterms is not known. However, one can show that at NNLO, HTLpt can 
be renormalizcd using only local counterterms for the vacuum, the Debye and fermion masses, and 
the coupling constant. The countcrterm for a s coincides with the perturbative value giving rise 
to the standard one-loop running. We do not list the counterterms, but present the full results 



elsewhere [llj. If the expansion in S could be carried out to all orders, the final result would be 
independent of the HTL parameters and m q . However, at any finite order in 5, the results 
depend on mo and m q . A prescription is then required to determine these parameters. We will 
discuss the prescription we use below. 



3. Thermodynamic potential 

In this section, we present the final results for the thermodynamic potential 51 at orders 5° 
(LO), S (NLO), and S 2 (NNLO). The LO and NLO results were first obtained in Refs. Q and 
they are listed here for completeness. At LO, the thermodynamic potential was calculated exactly, 
while at NLO and NNLO the resulting expressions for the diagrams arc too complicated. To make 
the calculations tractable, the thermodynamic potential is therefore evaluated approximately by 
expanding them in powers of mo/T and m q /T which assumes that these ratios are 0(g s ). This 
implies that the thermodynamic potential is evaluated in a double expansion in g s , mu/T, and 
rriq/T, and we have kept terms that contribute naively through order Due to the magnetic mass 
problem [l2| . HTLpt suffers from the same infrared divergences as ordinary perturbation theory 
and <7g is the highest order computable using only perturbative methods. 

The complete expression for the leading order thermodynamic potential is given by [H 

_ , , 7d F 15^ 2 on d F 2 3 45 { p. 7 



1 + ~ - -mjj - 30-^ + 30m^ + — log^ - - + lE + — 



-Fidcai 4 d A 2 u d A q u 4 V 2 2 3 

d A 



eo^V-e)^, (5) 
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where J-idcai = — — l)7r 2 T 4 /45 is the free energy of an ideal gas of noninteracting 
je is the Euler-Mascheroni constant. Moreover, we have introduced the dimensionless 
A = /x/27rT, mo = mo/^T, and m q = m q /2irT. 
The NLO thermodynamic potential reads [3] 
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Finally, our new result for the NNLO thermodynamic potential for QCD is 
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where is the Riemann zeta-function. 

As pointed out earlier, the HTL mass parameters are completely arbitrary and we need a 
prescription for them in order to complete a calculation. The variational mass prescription unfortu- 
nately gives rise to a complex Dcbye mass and m q = at NNLO. One strategy is therefore to throw 
away the imaginary part of the thermodynamic potential to obtain thermodynamic functions that 
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Figure 1: Comparison of LO, NLO, and NNLO predictions for the scaled pressure for Nf = 2 + 1 (left panel) and 
Nf = 2 + 1 + 1 (right panel) with lattice data from Cheng et al. 0] and Borsanyi et al. Q ■ We use N c = 3, three-loop 
running for a a . ft = 2+T, and A^r = 344 MeV. Shaded band shows the result of varying the renormalization scale 
fj, by a factor of two around fi = 2nT for the NNLO result. See main text for details. 



are real valued [131 14 1. Here we use another strategy explored in Rcfs. [H Il4j that is inspired by 



dimensional reduction: We equate the Debye mass with the mass parameter of three-dimensional 
electric QCD (EQCD) @, i. e. raj) =mg. In Ref. @, it was calculated to NLO giving 

47ra Srn2 f i c\a s ^5 11 11 n _ fi\ caSf&s 4 ] _ n 7 



7 fi s 2 F a s (l 4 2 2 fi\ 3 s 2F a s \ 

+ 6 l0g 2 + — lr3 l0g2 -3 7£ -3 l0g 2j'2 — ) " (8) 

This mass can be interpreted as the contribution to the Debye mass from the hard scale T and is 
well defined and gauge invariant order-by-order in perturbation theory. However, beyond NLO, it 
will also depend on factorization scale that separates the hard scale and the soft scale gT. For the 
quark mass, here we choose m q = 0. 

The final NNLO results are very insensitive to whether one chooses a perturbative mass pre- 
scription for m q or m q = 0; however, convergence is improved with the choice m q = 0. A detailed 
presentation of the full calculation of the NNLO thermodynamic potential and the dependence of 



our final results on the mass prescriptions for mu and m q is forthcoming in a longer paper 11 1 . 



4. Results 

In Fig. [1] we show the normalized pressure for iV c = 3 and Nf = 2 + 1 (left panel) , and N c = 3 
and Nf = 2+1 + 1 (right panel) as a function of T. The results at LO, NLO, and NNLO use 
the BN mass given by Eq. (j8|) as well as m q — 0. For the strong coupling constant a s , we used 
three-loop running [l| with = 344 MeV which for N f = 3 gives a s (5 GeV) = 0.2034 [11]. 
The central line is evaluated with the renormalization scale \i = lixT which is the value one expects 
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Figure 2: Comparison of NNLO predictions for the scaled trace anomaly with Nt = 2+1 (left panel) and Nf = 2+1+1 
fermions (right panel) lattice data from Cheng et al. 01 and Borsanyi et al. Q. We use N c = 3, three-loop running 
for a s , fj, = 2-kT, and A^jg- = 344 MeV. Shaded band shows the result of varying the renormalization scale fj, by a 
factor of two around fi = 2irT. See main text for details. 



from effective field theory calculations 0, and the band represents a variation of fi by a factor 
of two around this scale. 

The lattice data from the Wuppertal-Budapcst collaboration uses the stout action and have 
been continuum extrapolated by averaging the trace anomaly measured using their two smallest 
lattice spacings corresponding to iV r = 8 and iV r = 10 3]0 Using standard lattice techniques, 
the continuum-extrapolated pressure is computed from an integral of the trace anomaly. The 
lattice data from the hotQCD collaboration are their N T = 8 results using both the asqtad and 
p4 actions Q. The hotQCD results have not been continuum extrapolated and the error bars 
correspond to only statistical errors and do not factor in the systematic error associated with the 
calculation which, for the pressure, is estimated by the hotQCD collaboration to be between 5 - 
10%. We note that there are hotQCD results for physical light quark masses [lH; however, these 
are available only for temperatures below 260 MeV and the results are very close to the results 
shown in the figures so we do not include them here. 

As can be seen from Fig. [1] the successive HTLpt approximations represent an improvement 
over the successive app roximations coming from a naive weak-coupling expansion; however, as in 
the pure-glue case 14(, the NNLO result represents a significant correction to the LO and NLO 
results. That being said the NNLO HTLpt result agrees quite well with the available lattice data 
down to temperatures on the order or 2 T c ~ 340 MeV for both N f = 3 (Fig. [1] left) and N f = 4 
(Fig. Q] right). Below these temperatures the successive approximations give large corrections with 
the correction from NLO to NNLO reaching 100% near T c . 

In Fig. [2J we show the NNLO approximation to the trace anomaly (interaction measure) nor- 



1 We note that the Wuppertal-Budapest group has published a few data points for the trace anomaly with N T = 12 
and within statistical error bars these are consistent with the published continuum extrapolated results. 
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malizcd to T 4 as a function of T for N c = 3 and Nf = 3 (left panel) and for N c = 3 and JV/ = 4 
(right panel). In the left panel we show data from both the Wuppertal-Budapest collaboration 
and the hotQCD collaboration taken from the same data sets displayed in Fig. Q] and described 
previously. In the case of the hotQCD results we note that the results for the trace anomaly using 
the p4 action show large lattice size affects at all temperatures shown and the asqtad results for 
the trace anomaly show large lattice size effects for T > 200 MeV. In the right panel we display a 
parameterization (solid blue curve) of the trace anomaly for Nf = 4 published by the collaboration 
@] since the individual data points were not published. In both the left and right panels we see 
very good agreement with the available lattice data down to temperatures on the order of T ~ 2 T c . 



5. Summary and outlook 

We have presented results for the LO, NLO, and NNLO thermodynamic functions for SU(iV c ) 
Yang-Mills theory with Nf fermions using HTLpt. We compared our predictions with lattice data 
for N c = 3 and Nf £ {3, 4} and found that HTLpt is consistent with available lattice data down 
to T ~ 2 T c for the pressure and the trace anomaly. This is in line with expectations since one is 
expanding about the trivial vacuum — and therefore neglects the approximate center symmetry 
Z(N C ). Close to the deconfinement transition, it is essential to incorporate this symmetry (l9| . 



Comparing our results with the NNLO results of pure Yang-Mills [14J , we find that including 
the quarks gives much better agreement with lattice data. This is not unexpected since fermions 
are "perturbativc" in the sense that they decouple in the dimensional-reduction step of effective 
field theory. 

As was the case with pure Yang-Mills we found that the variational solution for the Debye 
mass mo is complex and we therefore chose instead to use the perturbative mass parameter from 
EQCD together with m q = 0. Whether the complexity of the variational Debye mass is due to the 
additional expansion in mu/T and m q /T is impossible to decide at this stage. We also found that 
there was a large correction going from NLO to NNLO. Unfortunately, due to the magnetic mass 
problem it is impossible to go to N 3 LO to see whether the problem persists without supplementing 
our calculation with input from three-dimensional lattice calculations. 

In closing, we emphasize that HTLpt provides a gauge invariant reorganization of perturba- 
tion theory for calculating static and dynamic quantities in thermal field theory. Given the good 
agreement with lattice data for thermodynamics, it would be interesting to apply HTLpt to the 
calculation of real-time quantities at temperatures that are relevant for LHC. 
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